Relaxation of Weakly Interacting Electrons in One Dimension 
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We consider the problem of relaxation in a one-dimensional system of interacting electrons. In the 
limit of weak interactions, we calculate the decay rate of a single-electron excitation, accounting for 
the nonlinear dispersion. The leading processes which determine the relaxation involve scattering of 
three particles. We elucidate how particular forms of Coulomb interaction, unscreened and screened, 
lead to different results for the decay rates and identify the dominant scattering processes responsible 
for relaxation of excitations of different energies. Interestingly, temperatures much smaller than the 
excitation energy strongly affect the rate. At higher temperatures the quasiparticle relaxes by 
exciting co-propagating electron-hole pairs, whereas at lowest temperatures the relaxation proceeds 
via excitations of both co-propagating and counter-propagating pairs. 

PACS numbers: 71.10.Pm 
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Low energy excitations of a three-dimensional interact- 
ing electron system are fermionic quasiparticles that in 
many respects resemble bare electrons [l|. A quasiparti- 
cle excitation of energy e has finite decay rate r _1 oc e 2 , 
where e is measured from the Fermi level. This fact is 
the foundation of the Fermi liquid theory and was con- 
firmed experimentally by measuring the broadening of 
the Lorenzian-shaped spectral function 

One-dimensional interacting fermions are convention- 
ally described within the exactly solvable Tomonaga- 
Luttinger model where particles are assumed to have 
a linear dispersion. This model can be diagonalized in 
terms of noninteracting bosonic excitations [3, |J] which 
have infinite lifetimes. This feature reveals an important 
limitation of the Tomonaga-Luttinger model, because in 
general an excited system is expected to relax to equi- 
librium. To study relaxation, one should therefore con- 
sider models that take into account the curvature of the 
spectrum. Recent experimental observation of different 
equilibration rates of hot electrons and holes in quan- 
tum wires [j| has confirmed the importance nonlinear 
dispersion of electrons. Study of interacting electrons 
with nonlinear spectrum is a subject of intense theoreti- 
cal interest [6H13j. The area of new physics beyond the 
Luttinger liquid formalism has been recently reviewed in 
Ref. 0- 

In this Letter we consider a system of spinless fermions 
and study the decay rate of a quasiparticle excitation 
placed above the Fermi level, see Fig. [TJ Since the com- 
plete study of the effects of nonlinear dispersion is very 
difficult, here we analyze the limit of weak interactions. 
In this case, the scattering processes can be classified 
by the number of colliding particles. Unlike in higher di- 
mensions, two-particle processes do not lead to relaxation 
due to the conservation laws of energy and momentum. 
Therefore, the leading mechanism which provides finite 
relaxation rate involves scattering of three particles [f|. 




FIG. 1. A single electron excitation of energy e and momen- 
tum pi relaxes via scattering off two other electrons. Filled 
and empty circles represent, respectively, incoming and out- 
going momenta. The one-dimensional topology of the Fermi 
surface determines the three cases (a), (b), and (c), which we 
refer to as the "2-1"-, "3-0"-, and "l-2"-type processes. While 
at zero temperature only the "2-1" processes are allowed, all 
three types contribute to the quasiparticle decay at nonzero 
temperatures. 



At zero temperature, the problem of relaxation was 
studied in Ref. Q, where the quasiparticle decay rate 
was found to behave as r _1 oc e 8 . Here we study the 
effect of temperature T, and find dramatic departures 
from the T = case. We consider the situation when 
the temperature is much smaller than the energy of the 
excitation. Denoting the momentum of the excitation by 
Pi, we find the expression for its decay rate 

- = £' W^ 3 'mn 3 (l - n v )(l - n 2 >){l - ny). (1) 



This equation accounts for scattering of the quasiparticle 
of momentum p\ and two others with momenta {p2,P3} 
into three outgoing states , P2 j } ; see Fig. [TJ In 
Eq. {TJ by rn = n Pi we denote the Fermi occupation num- 
bers, while £' indicates summation over distinct states. 
The scattering rate 3 is determined by the Fermi 
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golden rule expression, and reads 



with the even periodic function f c {&) — fc(S + t/3) 
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Here the delta function imposes conservation of the to- 
tal energy, defined as E — s Pl + e P2 + e P3 and similarly 
for the outgoing momenta. The three particle scattering 
amplitude is defined as a vacuum expectation value 
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It is the central object which determines the relaxation 
rate (|TJ) . The unperturbed Hamiltonian Hq and the per- 
turbation V are taken in the form 

H = ^2e p ala p , V = — ^ V q al i+q al 2 _ g a P2 a Pl) 
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where a and are the fermionic operators, e p = p 2 /2m, 
and L is the system size. The two-body interaction enters 
the Hamiltonian via its Fourier transform V q . 

The scattering amplitude ((3]) is very sensitive to the 
form of the two-body interaction V(x). It vanishes for 
V(x) oc 5(x), i.e., V q = const, corresponding to the con- 
tact interaction between spinless fermions. Nullification 
of the amplitude in this case arises because Pauli prin- 
ciple prevents two electrons from sharing the same po- 
sition. The amplitude© also vanishes for the Cheon- 
Shigehara model 3, 15j|7 defined by V{x) oc S"(x), i.e., 



V q oc q 2 . The latter belongs to the class of the so-called 
integrable models [lj| for which there is no relaxation 
and therefore t _1 = 0. 

In quantum wires, the interaction between electrons 
is of a longer range. As an example of the most prac- 
tical use, we consider the Coulomb interaction defined 
by V c (x) — e 2 /\x\, which has the Fourier transform 
2e 2 ln(l/|p|). Interestingly, this logarithmic form gives 
a vanishing three-particle amplitude (j3|), although the 
model describing fermions interacting via Coulomb in- 
teraction is not expected to be integrable. It is therefore 
important to cut off the short distance singularity. This 
is done by accounting for the finite width of the wire w. 
One then obtains 
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keeping the first two leading-order terms. Here and in 
the following we neglect numerical factors under the log- 
arithm because we consider small momenta \p\ <C fi/w. 

Calculation of the three particle amplitudes is rather 
tedious [ljj]. For unscreened Coulomb interaction (@}, 



one obtains the leading order result on the mass shell 
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Instead of using the incoming momenta {pj } and the out- 
going ones {p'j}, for convenience here we have introduced 
Jacobi coordinates P,p r ,ip, defined as 

p i =P/3+ v / 273prCos(^-2 7 ri/3) ) j = 1,2,3, (7) 

and similarly for the outgoing momenta. Here P — p\ + 
P2 + P3 has a simple meaning of the total momentum of 
the three colliding particles. The momentum p r is given 
by 



Pr 



Pi - P2 



) 2 + (P2- P3 ) 2 + (P3-Pi) 2 /V3, (8) 



and measures the typical separation between the mo- 
menta. We note that the conservation laws impose 

Pr =Pr- 

We consider scattering at low temperatures, when all 
scattering particles should be in the vicinity of the two 
Fermi points. This enables us to classify particles as 
being right- or left-moving. Throughout this paper we 
study the decay of a right-moving excitation. The two 
other incoming particles can be classified in one of three 
ways: (i) one particle has positive momentum and the 
other negative, (ii) both have positive momenta, and 
(iii) both have negative momenta. These three config- 
urations, respectively, have the values of the total mo- 
mentum near pp, 3pf, and ~pf, where pp is the Fermi 
momentum. Therefore, the incoming and the outgo- 
ing states must be in the same momentum configura- 
tion. Thus, we can distinguish between three different 
cases, which we call "2-1", "3-0", and "1-2" processes, 
see Fig. [1] We choose notations where the momenta pj 
and p'j are always on the same branch of the Fermi sur- 
face, while the scattering amplitude takes into account 
all possible exchange processes. 

At zero temperature, only processes of "2-1" type lead 
to relaxation [Tj. On the contrary, at nonzero tempera- 
tures all three processes have nonzero rates. In the follow- 
ing we will calculate the rates and identify the dominant 
processes as the temperature increases. 

We start our analysis by considering the "2-1" scatter- 
ing process, Fig.QJa). In order to understand the energy 
dependence of the decay rate, let us consider the mo- 
mentum change of the left-moving particle-hole pair. It 
is easily obtained from the conservation laws and reads 
P2 ~ P2 = (Pi ~ Pi) (pi - P3)/(P 3 - Pa). Using the esti- 
mate p% — P2 ~ %Pf, we can express the conservation 
laws contained in the scattering rate ([2]) as 
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where the momentum conservation comes from the am- 
plitude ©. We can now employ the last equation to per- 
form the summations over p' 2 and p' 3 in the expression for 
the rate (fTJ). The summation over the remaining three 
momenta P2,P3, and p' x then determines the rate. From 
Eq. ^ we conclude that the typical energy of the left 
moving particle-hole pair is of the order of e 2 /ep, where 
e ~ vp(pi —pf), and vf is the Fermi velocity. Therefore, 
the integration over p2 is restricted to an energy range of 
that width, while for both p% and p\ that range is of order 
e. As a result, at T = we find the phase space volume 
available for scattering to be proportional to e 2 (e 2 /ep), 
where ep denotes the Fermi energy Q. Since the ampli- 
tude ([5]) depends on momenta only logarithmically, we 
infer the scattering rate r _1 oc e 4 . This result still ap- 
plies at very low temperatures T <C e 2 /ep, since then the 
thermal smearing of the occupation numbers is not signif- 
icant. However, a new behavior of the decay rate emerges 
in the range of temperatures e 2 /ep -C T -C e, because 
the occupation numbers ri2,n2> in Eq. ([I]) of the left- 
moving pair become thermally smeared. As a result, the 
integration over momentum p2 covers the energy range of 
order T, and the phase space volume is proportional to 
e 2 T |14|. After careful calculation, using the amplitude 



calculation reveals the final result for the unscreened case 



([5]) and (O for the unscreened Coulomb interaction, one 
finds 
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where the numerical prefactors are c\ = (15 — 7r 2 )/10247r 3 
and C2 = 3(12 — 7r 2 )/647r 3 . The logarithm of the expres- 
sion ([TUl) originates from the first term in the amplitude 
(0, using p r — pp. 

At zero temperature the "3-0" processes, Fig. QJb), 
are not allowed by the conservation laws 0. However, 
they do contribute to the decay rate at finite temper- 
atures. To evaluate the rate one can employ a simi- 
lar strategy to the one used above for the "2-1" pro- 
cesses. The momentum change of the initial excitation is 
Pi - Pi = (p3 ~ P3)(P3 " ftO/bi ~ P2), which enables us 
to express the conservation laws as 
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We can now use Eq. (1111) to perform the summation over 
the momenta p[ and P2 in Eq. (TTJ), which gives rise to 
a factor of 1/e in the rate. The remaining summation 
over P21P31P3 is over the typical range of momenta of 
the order T/vp and delivers the factor T 3 . The detailed 
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where C3 f=s 1.13 [17]]. The logarithmic prefactor arises 
from the amplitude ([3]), using p r ~ pi — pp. 

The processes of "1-2" type are similar to the "3-0" 
ones. In order to evaluate their decay rate, one can use 
the expression for the momentum change p\ — p[ of the 
initial excitation of the same form as for the "3-0" case. 
However, in the present situation the denominator p\ —p'2 
should be replaced by 2pp, rather than by p\— pf- There- 
fore, the final result for the rate should be the same as in 
Eq. (fT2"j) . provided one replaces e by iep, which we con- 
firmed by a careful calculation. Since e <C ep, the con- 
tribution of the"l-2" processes is always subdominant. 

We now turn to the case of screened Coulomb inter- 
action. We model the two-body potential as V sc (x) — 
e 2 /\x\ — e 2 j\Jx 2 + Ad 2 , where d is the distance between 
the wire and a conducting plane representing nearby 
gates. At |p| <C h/d <C h/w, its Fourier transform is 



K sc = 2e 2 In f - ) - 2e 2 



The first term in the right hand side of Eq. ([131) corre- 
sponds to the contact interaction V oc S(x), which does 
not affect spinless fermions. Therefore, the three-particle 
amplitude is determined by the second term of Eq. (|13[) 
and reads [l?} 
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Compared to the unscreened interaction (j4j) , the screened 
interaction (|13j) has two additional powers of momentum. 
This is reflected in the corresponding amplitude (|14|) , up 
to the logarithmic terms. As a result, the decay rates 
will have four additional powers of energy with respect 
to the unscreened case. For the "2-1" processes, the char- 
acteristic energy change is e, which determines the decay 
rate at zero temperature t -1 oc e i 
calculation one obtains 



Q. After a careful 
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FIG. 2. Leading behavior of the decay rate of a quasiparticle 
of energy e for (a) unscreened Coulomb interaction Q, and 
for (b) screened Coulomb interaction (|13[) . 



where c 4 = 15/32768tt 3 and c 5 = 55/9216tt 3 . For the 
"3-0" processes, the typical momentum change is T/vp 
and therefore one obtains four additional powers of tem- 
perature compared to Eq. (fT2l , 



ergy scale e* ~ (T 2 ^) 1 / 3 . Quasiparticles of energies 
lower than e* decay by exciting co-propagating particle- 
hole pairs, while quasiparticles of energies larger than 
e* decay by exciting both co-propagating and counter- 
propagating pairs. The same general picture applies in 
the case of the screened Coulomb interaction (ITBl . but 
the crossover energy scale e* is of order (T 6 ^) 1 / 7 . 

To summarize, motivated by recent experiment [E| we 
have calculated the decay rate of quasiparticles in weakly 
interacting one-dimensional electron systems. The dom- 
inant mechanism of quasiparticle decay involves three 
electrons and is illustrated in Figs.QJa) andQJb). The de- 
cay rate shows nontrivial temperature dependence even 
at T < e, see Fig. [2 
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where eg ~ 22.1. The decay rates given by Eqs. (|T0|) and 
(fT2")l for the Coulomb interaction as well as Eqs (fTS)) and 
(|17p for the screened Coulomb interaction are the main 
results of this Letter. Now we comment on their physical 
meaning. 

The spectral function of a system of interacting elec- 
trons described by the Tomonaga-Luttinger model dis- 
plays a power-law edge singularity on the mass shell 
m\l4. This diver gence is a signature of the infinite life- 
time of excitations. However, once one accounts for the 
curvature of the spectrum, the divergence disappears and 
the spectral function becomes broadened Q. Therefore, 
the quasiparticles on the mass shell are subject to decay. 
The above-calculated decay rate describes broadening of 
the spectral function in the vicinity of the particle mass 
shell. Our result (1T6|) taken at zero temperature is con- 
sistent with Ref. [3|. It is worth noting, however, that 
unlike our paper, where for the screened interaction all 
the relevant momentum scales are assumed to be small 
compared with h/d, Ref. 0] assumes pp ^> h/d, i.e., the 
Fourier components of interaction potential for momenta 
of the order of the Fermi momentum were neglected. 

The leading behavior of the decay rate is summa- 
rized in Fig. [5J Excitations of energies much larger than 
(Tap) 1 / 2 decay with a temperature independent rate. 
For the unscreened interaction (|4]), we infer a new en- 
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